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ABSTRACT 

We study the extent and covering fraction of cool baryons around galaxies of different luminosity and mass, 
based on a survey of Mg II AA2796,2803 absorption features near known galaxies. The initial sample consists 
of 13 galaxy and absorber pairs and 10 galaxies that do not produce Mg II absorption lines to within sensitive 
upper limits. The redshifts of the galaxy and absorber pairs range from z = 0.2067 to 0.892 with a median 
of z = 0.3818. We find that galaxies at larger impact parameters produce on average weaker Mgll absorbers. 
This anti-correlation is substantially improved when accounting for the intrinsic luminosities of individual 
galaxies. In addition, there exists a distinct boundary at p = /? gas , beyond which no Mgll absorbers are found. 
A maximum likelihood analysis shows that the observations are best described by an isothermal density profile 
and a scaling relation fl gas = 91 x (L B /L B J (035±0X>5) h~ l kpc (or 69 h~ l kpc at W(2796) = 0.3 A) with a mean 
covering factor of (n) = 80- 86 %. Together with the scaling relation between halo mass and galaxy luminosity 
inferred from halo occupation studies, this scaling of R gsis indicates that gas radius is a fixed fraction of the dark 
matter halo radius. We compare our results with previous studies and discuss the implications of our analysis 
for constraining the baryon content of galactic halos and for discriminating between competing scenarios for 
understanding the nature of the extended gas. 

Subject headings: cosmology: observations — intergalactic medium — quasars: absorption lines 



1, INTRODUCTION 

The Mg II AA2796,2803 doublets are among the absorption 
features commonly seen in the spectra of distant quasars that 
are produced by intervening gaseous clouds along the quasar 
lines of sight. Their rest-frame absorption equivalent width 
W(2796), ranging from W(2796) < 0.3 A to W(2796) > 2 
A (e.g. Steidel & Sargent 1990; Nestor et al. 2005; Prochter 
et al. 2006a), is found to represent the underlying gas kine- 
matics (e.g. Petitjean & Bergeron 1990; Churchill et al. 2000). 
Based on comparisons of the abundance ratios of various as- 
sociated ions, these Mgll absorption transitions are under- 
stood to arise primarily in photo-ionized gas of temperature 
T ~ 10 4 K (Bergeron & Stasmska 1986; Hamann 1997) and 
neutral hydrogen column density A^(HI) = 10 18 - 10 22 cm" 2 
(Churchill et al. 2000; Rao et al. 2006). In addition, surveys 
for galaxies near the observed Mg II absorbers have often un- 
covered luminous galaxies at projected distances (p < 50 h~ l 
kpc) and velocity separations (Av < 250 kms -1 ) from the 
absorbers (Bergeron 1986; Lanzetta & Bowen 1990, 1992; 
Steidel et al. 1994; Zibetti et al. 2005; Nestor et al. 2007; 
Kacprzak et al. 2007). The relatively large associated A^(HI) 
and the small separation between Mg II absorbers and galax- 
ies along common lines of sight indicate that Mg II absorbers 
may offer a sensitive probe of photo-ionized halo gas around 
galaxies at redshifts from z = 0.3 to z = 2.3 in the optical spec- 
tral window. 

Understanding the physical origin of these Mg II absorbers 
bears significantly on all efforts to apply their known statis- 
tical properties for constraining the baryon content of dark 
matter halos on different mass scales (Tinker & Chen 2008). 
Various theoretical models have been developed in the past 
that describe the nature of the absorbers in the context of 
gas accretion, including ram pressure stripped gas from ac- 
creted satellites (Wang 1993), gravitationally bound cold gas 
in halo substructures (e.g., Sternberg et al. 2002), and con- 



densed cold clouds in a thermally unstable hot halo (Mo & 
Miralda-Escude 1996; Mailer & Bullock 2004; Chelouche 
et al. 2007). 

Insights to the origin of Mgll absorbers can be obtained 
from their clustering amplitude on large scales (> 1 h~ l Mpc). 
The clustering of the absorbers is a consequence of dark 
matter halos in which they are found. Measurements of 
the galaxy and absorber cross-correlation function offer a 
means of quantifying the mean halo mass of the absorbers. 
High-mass halos are expected to be highly clustered, while 
low mass halos on average have weaker clustering strength. 
The cross-correlation function of Mg II absorbers and lumi- 
nous red galaxies (LRGs) have been measured at (z) = 0.6 
by Bouche et al. (2006), who found that Mg II absorbers of 
W(2796) > 0.3 A cluster strongly with LRGs and that weaker 
absorbers with W(2796) = 0.3- 1.15 A on average arise in 
dark matter halos that are 10 times more massive than those 
of stronger absorbers with W(2796) = 2-2.85 A. 

Both the observed clustering amplitude and the inverse cor- 
relation between the mean halo mass and absorber strength 
are difficult to interpret. For example, the frequency distri- 
bution function of Mg II absorbers show that there are on av- 
erage 10 times more W(2796) = 1 A absorbers than those of 
W(2796) = 2 A (e.g. Nestor et al. 2005; Prochter et al. 2006a) 
which, when combined with the clustering measurements, im- 
plies that the majority of absorbers reside in massive, highly 
biased halos. At the same time, we expect that observed Mg + 
ions originate primarily in photo-ionized gas of temperature 
T ~ 10 4 K and the halo gas in massive dark matter halos be- 
comes too hot for abundant Mg + to survive. 

In Tinker & Chen (2008; hereafter TC08), we developed 
a new technique that adopts the halo occupation framework 
for understanding the origin of QSO absorption-line systems. 
Specifically, the technique adopts a model density profile for 
the Mg + ions in individual dark matter halos. The "cold" 
baryon content of individual dark matter halos, as probed by 
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the presence of Mg + , is then constrained through matching 
the space density and clustering amplitude of dark matter ha- 
los with the observed frequency distribution function of Mg II 
absorbers and their clustering amplitude. Our model allows 
the possibility that a predominant fraction of the gas in mas- 
sive halos is shock heated to the virial temperature of the halo 
and becomes too hot to host abundant Mg II absorbing clouds. 
Within the hot halo, we further allow the possibility that some 
dense, cold clouds may penetrate through as seen in high res- 
olution cosmological simulations (Kravtsov 2003). The re- 
sult of this halo occupation analysis is the probability function 
P(W\Mh) that characterizes, for each dark matter halo of Mh, 
the total probability of finding a Mg II absorber of equivalent 
width W. 

The halo occupation analysis shows that observations, in- 
cluding both the Mg II absorber frequency distribution func- 
tion and their clustering amplitude, demand a rapid transition 
in the halo gas content at M"' 1 ' ~ 10 11 5 h~ l M Q . Below M c h "', 
halos contain predominantly cold gas and therefore contribute 
significantly to the observed Mgll statistics. Beyond Mf u , 
the cold gas fraction is substantially reduced and presumably 
the halo gas becomes too hot to maintain a large contribution 
to Mg II absorbers. In order to reproduce the observed overall 
strong clustering of the absorbers and the inverse correlation 
between W(2796) and halo mass Mh, roughly 5% of the gas in 
halos up to 10 14 h~ x M is required to be cold. It is understood 
under our model that the W(2796) vs. Mh inverse correlation 
arises mainly as a result of an elevated clustering amplitude 
of W (2796) < 1 A absorbers, rather than a suppressed clus- 
tering strength of "W '(2796) > 2 A absorbers. The clustering 
amplitude of weak absorbers is elevated because of the pres- 
ence of cold streams in massive hot halos. The amount of 
cold gas in clusters is small, therefore these halos can only 
contribute to weak absorbers. 

The initial results of our halo occupation analysis demon- 
strate that combining known statistics of dark matter halos 
with a simple model for their gas content can already repro- 
duce the statistical properties known for Mgll absorbers. It 
provides a simple prescription for populating baryons in dark 
matter halos that can be compared directly to results from nu- 
merical simulations and offer insights for understanding the 
physics of gas accretion in halos of all mass scales. However, 
some uncertainties remain in the halo occupation analysis. 

First, a key parameter that constrains the distribution of 
Mg + ions in individual dark matter halos is their incidence 
rate, K(Mf,), per halo. It specifies the total probability of de- 
tecting an absorber in a halo of mass M/,. In the halo occu- 
pation analysis of TC08, we find that dark matter halos of 
Mi, = io 11 ' 5 " 12 ' 5 h~ l M© on average have 100% incidence rate 
of Mg + ions, (n(M h )) = 1, at R gns < /?20o/3 for Mg + ions and 
(n{Mh)) = at larger radii. 1 In halos of lower masses, (k(M/,)) 
declines sharply. 

While a 100% incidence rate shows that every dark mat- 
ter halo hosts a uniform gaseous halo of Mg + ions with size 
/? gas , the low (/t(Mft)) is more difficult to interpret because 
(k(Mi,)) represents a mean value averaged over all halos of 
mass M/,. A low (k(M/,)) may be a result of only a small frac- 
tion of dark matter halos containing extended distributions of 
Mg + ions or a result of a small covering factor of Mg + in all 

1 The halo size #200 corresponds to the radius, within which the enclosed 
mean density is 200 X p m and p m is the background density. It is motivated 
by estimates of the virial radius both from the spherical collapse model and 
from N-body simulations that predict ss 180 X p m . 



halos. In addition, we have assumed in our initial analysis 
that the gaseous extent is related to halo mass according to 
flgas = [M h /(W n hr l M Q )f, and fl gas » = 50 h~ l physical 
kpc and (3 = 1/3. This scaling relation follows the theoretical 
expectation between virial radius and halo mass, but a higher 
^gas would naturally lower k for a fixed frequency distribution 
function of Mg II absorbers. 

Second, the physical origin of the cold clouds probed by the 
Mgll absorption transitions is ambiguous. While the mass 
scale at which halo gas is found to experience a transition 
from a cold-mode dominated state to a hot-mode dominated 
state agrees well with the expectation of theoretical models 
established to characterize the gas accretion history in dark 
matter halos (e.g. Birnboim & Dekel 2003; Keres et al. 2005; 
Dekel & Birnboim 2006; Birnboim et al. 2007), the observed 
cold-hot transition can also be interpreted as a declining for- 
mation efficiency of cool clouds in a hot halo with increas- 
ing halo mass (e.g. Mo & Miralda-Escude 1996; Mailer & 
Bullock 2004). Both of these scenarios predict that halo gas 
is primarily cold at M h < 10 11 5 /z _1 M . On the other hand, 
Bouche et al. (2006) have argued that the inverse correlation 
between W(2796) and clustering amplitude is suggestive of 
the absorbers arising in non-virialized gas flows, such as star- 
burst winds. 

While dense clumps in starburst driven outflows are ex- 
pected to contribute to some fraction of the observed Mg II 
absorbers, an important quantity to specify is the significance 
of this fraction. If a large fraction of Mg II absorbers origi- 
nate in cold dense clumps in starburst winds, then one must 
also account for the more complex star formation physics in 
efforts to apply known Mg II statistics for constraining the gas 
content of dark matter halos. We defer to § 5.3 for a more 
detailed discussion of the caveats in interpreting the known 
properties of Mg II absorbers as evidence to support their ori- 
gin in starburst winds. 

Direct constraints on the covering fraction of Mg II absorb- 
ing gas versus halo mass not only help to break the degen- 
eracy between the fraction of dark matter halos containing 
cold baryons and the extent of cold gas in individual halos, in 
principle they also serve to discriminate between the starburst 
wind scenario and gas accretion. Under the starburst wind 
scenario, superwind driven outflows are expected to proceed 
over a finite angular span along the minor axis (e.g. Heckman 
etal. 1990; Veilleux et al. 2005), and result in only partial 
covering of the halos. In addition, because at a given epoch 
only some fraction of galaxies are found at starburst or post- 
starburst stages, we would expect a fraction of dark matter ha- 
los to contain extended Mg II absorbing gas due to outflows. 

To obtain empirical constraints on the extent and cover- 
ing fraction of Mg + ions around galaxies of different lumi- 
nosity and mass, we have initiated a survey of Mgll ab- 
sorbers around known galaxies at small projected distances 
(p < 100 h~ l kpc) to a background QSO. Accounting the pres- 
ence or absence of Mg II absorbers for an unbiased sample of 
galaxies at different impact parameters yields a statistical esti- 
mate of the gas covering fraction around field galaxies. A de- 
tailed comparison between the absorber strength and galaxy 
properties further constrains the density profile and extent of 
the gas. The primary objectives of our study are (1) to im- 
prove the uncertainties of our halo occupation analysis by 
obtaining empirical measurements of R g . ds and (3, and (2) to 
examine whether starburst driven outflows are the predomi- 
nant mechanism for producing the observed Mg II absorbers 
by obtaining an empirical constraint of the gas covering factor 
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k. Here we present the initial results of our study. 

This paper is organized as follows. In Section 2, we de- 
scribe the design of our experiment for constraining the ex- 
tent and covering fraction of Mg + ions. In Section 3, we de- 
scribe the imaging and spectroscopic data available for our 
experiment. In Section 4, we examine the correlation be- 
tween Mg II absorption strength and galaxy properties, such 
as the impact parameter and luminosity. We compare our 
results with previous studies and discuss the implications of 
our analysis in Section 5. We adopt a ACDM cosmology, 
Om = 0.3 and Q\ = 0.7, with a dimensionless Hubble constant 
h = //(>/( 100 km s" 1 Mpc" 1 ) throughout the paper. 

2. EXPERIMENT DESIGN 

Accurate estimates of the extent and covering fraction of 
Mg + ions around galaxies require an unbiased survey of Mg II 
absorbers in the vicinity of known galaxies. Nearly all studies 
published in the literature to date focus on surveys of galax- 
ies that correspond to known Mg II absorbers (e.g. Bergeron 
1986; Lanzetta & Bowen 1990, 1992; Steidel et al. 1994; Zi- 
betti et al. 2007; Kacprzak et al. 2007), with the exception 
of an on-going study described in Tripp & Bowen (2005). 
Searching for galaxies that give rise to known Mg II absorbers 
bears a large uncertainty due to contamination from galaxies 
that are correlated with the true absorbing galaxies. In addi- 
tion, these surveys based on known presence of Mg + ions do 
not yield constraints on the incidence of these ions in galactic 
halos. 

We are conducting a survey of Mg II absorbers at the lo- 
cations of known galaxies close to the line of sight toward a 
background QSO. We identify galaxies at impact parameters 
p < ro. The value of ro is selected to allow identifications 
of galaxies at impact parameters beyond the size of extended 
Mg II halos, in order to constrain the extent and covering frac- 
tion of the gas. At the same time, we require ro to be small 
enough in order to minimize the contamination of galaxies 
that do not give rise to the absorbers but are merely correlated 
with the absorbing galaxies through clustering. 

To estimate the contamination rate of correlated galaxies 
as a function of projected distance, we adopt the distribution 
of dark matter halos identified in cosmological simulations as 
a guide. We use the z = 0.5 output of a 400 h~ l co-moving 
Mpc simulation described in Tinker et al. (2007). First, we 
identify for each random line of sight dark matter halos that 
have Mh > 10 11 /z _1 M Q and are located at p < /?20o/3 from 
the sightline. From TC08, we consider these halos likely to 
host a Mg II absorber along the line of sight. This is also con- 
sistent with the results from the halo occupation distribution 
models of galaxies found in DEEP2 and Sloan Digital Sky 
Survey from Zheng, Coil, & Zehavi (2007; hereafter ZCZ07). 
These authors found that the minimum mass of a dark matter 
halo that on average hosts a galaxy of > 0.1L, is w 2 x 10 11 
h~ l M©. Next, for each potential Mg II absorber, we search for 
neighboring halos that are at \5v\ < 250 kms" 1 and projected 
distance less than p from the absorber and may therefore con- 
taminate the search of the true absorber. Finally, we repeat 
this exercise over 10,000 random lines of sight and calculate 
the frequency of finding contaminating halos as a function of 
P- 

The result of our exercise is presented in Figure 1, where 
we show the probability of contamination due to correlated 
galaxies versus impact parameter. The errorbars represent the 
68% scatter of the mean from 16 jackknife subsamples. Fig- 
ure 1 shows that at p > 100 h~ l physical kpc there is a > 30% 
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FIG. 1. — The probability of finding more than one galaxy with <5v < 250 
kms~' and at projected distances less than p from an absorber. The errorbars 
represent the 68% scatter that was estimated using a jackknife resampling 
technique. 



chance that one would mis-identify a correlated galaxy as the 
true absorber. Combining the contamination rate estimates 
with previous search reports that indicate a mean gaseous ex- 
tent of ~ 50 h~ l physical kpc for luminous galaxies (Bergeron 
1986; Lanzetta & Bowen 1990, 1992; Steidel et al. 1994), we 
carry out our study based on a sample of galaxies that are 
within ro = 100 h~ x kpc from a background QSO. 

3. DATA 

To study the extent and covering fraction of Mg + ions in 
galactic halos, it requires survey data of both galaxies and 
absorbers along common lines of sight toward background 
quasars. We have selected seven QSO fields for which galax- 
ies data are already available, and echelle spectra of the back- 
ground QSOs have also been obtained either from our own 
observations or through the ESO science data archive. We de- 
scribe available QSO spectra and galaxy data separately in the 
following sections. 

3.1. Echelle Spectra of QSOs 

Echelle spectroscopic observations of the QSOs, 
PKS 0122-0021 (zqso = 1.0765), HE 0226-4110 
( ZQSO = 0.495), PKS 0349-1438 (z QS0 = 0.6163), and 
PKS 0454-2203 (z QS0 = 0.5335) were obtained using the 
MIKE echelle spectrograph (Bernstein et al. 2003) on the 
Magellan Clay telescope in 2007 October. The spectrograph 
contains a blue camera and a red camera, allowing a full 
wavelength coverage from near-ultraviolet 3350 A through 
near-infrared 9400 A. Depending on the brightness of the 
QSO, the observations were carried out in a sequence of one 
to three exposures of duration 600 s to 900 s each. The mean 
seeing condition over the period of integration was 0.7". We 
used a 1" slit and 2x2 binning during readout, yielding a 
spectral resolution of FWHM w 13 kms" 1 at wavelength 
A = 4500 A and w 16 kms" 1 at A = 8000 A. The data 
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were processed and reduced using the MIKE data reduction 
software developed by Buries, Prochaska, & Bernstein 2 . 
Wavelengths were calibrated to a ThAr frame obtained 
immediately after each exposure and subsequently corrected 
to vacuum and heliocentric wavelengths. Flux calibration 
was performed using a sensitivity function derived from 
observations of the flux standard EG131. The flux calibrated 
spectra have a mean S/N > 25 per resolution element across 
the entire spectral range. 

Echelle spectra of the QSOs, PKS1354+1933 (z QS0 = 
0.7200), PKS1424-1150 (z QS o = 0.8060), and PKS1622+23 
(zqso = 0.927) were obtained using the Ultraviolet and Vi- 
sual Echelle Spectrograph (UVES; D'Odorico et al. 2000) on 
the VLT/UT2 telescope and were retrieved from the ESO 
data archive. Observations of PKS1354+1933 (program ID 
076.A-0860) and PKS 1424-1 150 (program ID 075.A-0841) 
were obtained using a 1" slit and the dichroics #1 B346+R580 
and #2 B437+R860 in a sequence of two to three exposures 
of duration 180s to 300 s each. The data were binned 2x2 
during readout. Observations of PKS 1622+23 (program ID 
069.A-0371) were obtained using a 1" slit and the dichroic 
#1 B346+R580 in a sequence of three exposures of duration 
4900 s each. The data were also binned 2x2 during read- 
out. The UVES spectra were processed and calibrated using 
the standard data reduction pipeline. The combined spectra of 
the three QSOs cover a spectral range from 3300 A to 6600 
A, with a pixel scale of 0.035 A per pixel and S/N > 10 per 
pixel at A > 3700 A. The high S/N and high spectral resolu- 
tion MIKE and UVES spectra together allow us to search for 
Mgll absorption features at redshifts 0.2 < z < 2.3 to a lim- 
iting absorption strength W(2796) = 0.01 A along the lines of 
sight toward the seven background QSOs. 



3.2. The Galaxy Sample 

For the purpose of the study described in § 2, we need to 
identify a representative sample of galaxies at p < 100 ft -1 
kpc from the background QSOs. Galaxies in six of the seven 
QSO fields listed in Table 1 have been observed and ana- 
lyzed by Chen et al. (1998, 2001a). Specifically, galaxies 
in the fields around PKS 0122-0021, PKS 0349-1438, and 
PKS 0454-2203, PKS1354+1933, and PKS1424-1150 were 
observed as part of the survey program to study the gaseous 
extent of galaxies based on the presence or absence of cor- 
responding Lyct and C IV absorption features in the spectra 
of the background QSOs (Lanzetta et al. 1995; Chen etal. 
1998, 2001a; Chen et al. 2001b). The authors carried out a 
redshift survey of galaxies brighter than R = 21.5 mag at an- 
gular distances 9 < 1.3'. The survey depth enabled findings 
of galaxies with luminosities greater than 0.01 -1.0 L, and 
impact parameters less than 100-350 h~ l kpc at the redshifts 
0.1 < z ga ] < 0.6. The galaxy sample is therefore ideally suited 
for our present study. 

Galaxies in the field of PKS 1622+2352 were first published 
in Steidel et al. (1997) as part of a program to determine the 
redshifts of galaxies found in deep Hubble Space Telescope 
(HST) images that are brighter than R « 24.5. Images of in- 
dividual galaxies were also analyzed by Chen et al. (1998). 
Galaxies in the field of HE 0226-41 10 were observed in an 
on-going program to uncover galaxies brighter than R = 23 

2 see http://web.mit.edu/~burles/www/MIKE/mike_cookbook.html. 



mag at z ga i < 0.5 for understanding the cross-correlation be- 
tween galaxies and O VI A A 1031, 1037 absorbers (Chen et al. 
2008 in preparation). The two galaxy samples satisfy the se- 
lection criterion of an unbiased galaxy survey in QSO fields 
that does not rely on known knowledge of the presence of 
an Mg II absorber. We therefore include the galaxy data in 
our sample. Images of the galaxies are presented in Figure 
2. The data were obtained using the Wide Field and Plane- 
tary Camera 2 (WFPC2) and the F702W filter (see also Chen 
et al. 2001a) on board HST, with the exception of the images 
of galaxies in the field around HE 0226-41 10 that were ob- 
tained using IMACS and the B-band filter on the Magellan 
Baade telescope. 

The sample of galaxies identified at p < 100 h~ l kpc from 
the lines of sight toward seven background QSOs is pre- 
sented in Table 2, where we list for each galaxy in columns 
(1) — (7) the field, Right Ascension and Declination offsets 
from the QSO Aa and AS, redshift z ga i, impact parameter 
p, apparent B-band magnitude, and absolute B-band mag- 
nitude M B -5\ogh? Measurement uncertainties in R and 
Mb -5 log h were typically 0.2 dex. Our initial galaxy sam- 
ple consists of 23 galaxies at redshifts from z ga i = 0.2067 to 
Zgai = 0.892 with a median of z ga i = 0.3818 4 . We have ex- 
cluded galaxies that are within |Av| < 3000 kms" 1 from the 
background QSOs, where most of the gas is expected to be 
highly ionized due to the QSO proximity effect (e.g. Bajt- 
lik et al. 1988). The rest-frame B-band magnitudes of the 
galaxies range from M#-5 log h = -15.9 to -20.3 with a me- 
dian of M B -5 log h = -19.0, corresponding to 0.5L Bt for 
Mb, -5 log h = -19.8 that characterizes the blue galaxy pop- 
ulation at z ga i ~ 0.4 (e.g. Faber et al. 2007). 

3.3. The Galaxy-Mg II Absorber Pair Sample 

To establish a galaxy-Mg II absorber pair sample for the 
subsequent analysis, we examine the echelle spectra of the 
background QSOs and search for the corresponding Mg II ab- 
sorption doublet at the locations of the galaxies presented in 
§ 3.2. Specifically, we accept absorption lines according to a 
3cr detection threshold criterion, which is appropriate because 
the measurements are performed at a small number of known 
galaxy redshifts. Next, we form galaxy and absorber pairs by 
requiring a velocity difference 5 A v < 250 kms -1 . Finally, we 
measure 3-er upper limits to the 2796-A absorption equivalent 
widths of galaxies that are not paired with corresponding ab- 
sorbers. The procedure identifies 13 galaxy and absorber pairs 
and 10 galaxies that do not produce Mgll absorption lines to 
within sensitive upper limits. Impact parameter separations of 
the galaxy and absorber pairs range from p= 16.3 to 100.1 h~ l 
physical kpc with a median of p = 54.8 h~ l kpc. The corre- 
sponding Mg II absorption strength of each of the 23 galaxies 
in our sample is listed in columns (8) — (9) of Table 2. 



3 All magnitudes presented in this paper are in the AB system. 

4 We note that 12 of the 23 galaxies in our sample are from the field around 
3C336. These 12 galaxies have spectroscopic redshifts evenly distributed 
from z = 0.261 to z = 0.892 with a median of (z) = 0.6, covering a comoving 
distance of as 1 .4 Gpc. Therefore these galaxies are not part of a coherent 
structure. There is no known galaxy clusters along the sightline toward 3C336 
and there is no physical mechanism that suggests that slight overdensities on 
scales of 0.1-1 Gpc would affect gas properties within individual halos. 

5 We note that the galaxy at z ga i = 0.4592 in the field around 1354+1933 
is R2 550 kms from the identified Mgll absorber. However, including the 
redshift measurement uncertainty of 300 km s~' reported in Ellingson et al. 
(1991), we consider this pair a match. 
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TABLE 1 

Summary of the Echelle Spectroscopic Observations of Distant QSOs 



Field RA(J2000) Dec(J2000) zq S0 V Instrument Exptime Spectral Coverage 



PKS 0122-0021 01:25:28.8 -00:05:56 1.0765 16.4 MIKE 1200 3350-9400 

HE0226-4110 02:28:15.2 -40:57:16 0.4950 15.2 MIKE 600 3350-9400 

PKS 0349-1438 03:51:28.5 -14:29:09 0.6163 16.2 MIKE 1800 3350-9400 

PKS 0454-2203 04:56:08.9 -21:59:09 0.5335 16.1 MIKE 2700 3350-9400 

PKS 1354+1933 13:57:04.4 +19:19:07 0.7200 16.0 UVES 600 3300-6800 

PKS 1424-1150 14:27:38.1 -12:03:50 0.8060 16.5 UVES 720 3300-6600 

3C336 16:24:39.1 +23:45:12 0.9274 17.5 UVES 14,700 3300-6600 
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FIG. 2. — Images of the 23 galaxies in our sample that are identified at p < 100 h' 1 kpc from the lines of sight toward a background QSO. With the exception of 
the two galaxies in the field of HE 0226-4 110, the galaxy images were obtained using WFPC2 and the F702W filter (see also Chen et al. 2001a). Images of the 
galaxies in the field of HE 0226-41 10 were obtained using IMACS and the R filter on the Magellan Baade telescope. The spatial extent of each image is roughly 
25 fr 1 physical kpc on a side at the redshift of the galaxy indicated in the upper right-corner. The orientation of each image is arbitrary. The absorption strength 
of the corresponding Mgll absorbers found in the spectra of the background QSO described in § 3.1 is indicated in the lower-right corner. For non-detections, 
we provide a 3-cr upper limit to the rest-frame absorption equivalent width. 



4. THE EXTENT AND COVERING FACTOR OF MG+ IONS IN 
GALACTIC HALOS 

In this section, we examine the correlation between galaxies 
and Mg II absorbers identified along common lines of sight, 
using the sample of 13 galaxy-Mg II absorber pairs and 10 
galaxies at p < 100 h~ x kpc that do not give rise to Mgll ab- 
sorption to a sensitive upper limit. The choice of j-q = 100 h~ l 
kpc in selecting galaxies at p < ro for our analysis is adopted, 



in order to minimize the contamination of galaxies that are 
correlated with the true Mgll absorbing galaxies through 
clustering and to maximize the fraction of physical galaxy and 
absorber pairs in our sample. In addition to the simulation re- 
sults presented in § 2, this criterion is further supported by the 
stronger clustering signal at small (p < 100 h~ l kpc) impact 
parameters of LRGs (Bouche et al. 2004). This strong signal 
is a natural consequence of cross-correlating objects inside 
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TABLE 2 

Galaxies and Absorption Systems 



Galaxies Absorption Systems 





Aa 


AS 




p 




M B 




W(2796) 


Field 


(arc sec) 


(arc sec) 


-iial 


(h~ l kpc) 


R a 


— 5 log h 




(Al 


V 1 1 


(2) 


(3) 


(4) 


(5) 


(6) 


(1) 


(8) 


(9) 


0122-0021 


-8.7 


-12.3 


0.3788 


54.8 


20.6 


-19.2 


0.3791 


0.05 ±0.01 


0226-4110 


-7.9 


-7.9 


0.2067 


26.5 


22.2 


-16.4 




<0.02 




11.6 


-12.5 


0.2674 


49.0 


20.4 


-18.4 


0.2678 


0.03 ±0.01 


0349-1438 


-9.4 


10.8 


0.3567 


50.1 


20.6 


-19.1 


0.3572 


0.21 ±0.01 




11.7 


-24.1 


0.3236 


88.0 


20.0 


-19.1 




<0.01 


0454-2203 


10.6 


5.7 


0.2784 


35.6 


21.1 


-18.0 




<0.01 




-1.1 


-18.0 


0.4847 


75.8 


20.1 


-20.3 


0.4834 


0.42 ±0.01 




0.3 


-19.7 


0.3818 


72.0 


20.8 


-19.1 




<0.01 


1354+1933 


1.2 


7.5 


0.4592 


31.0 


21.0 


-19.3 


0.4563 


0.85 ±0.01 




-21.6 


-12.2 


0.4406 


98.8 


20.9 


-19.1 




<0.01 


1424-1150 


-0.2 


17.6 


0.3404 


59.8 


20.5 


-19.1 


0.3415 


0.12 ±0.04 


3C336 


3.0 


-0.0 


0.8920 


16.3 


23.3 


-19.0 


0.8909 


1.55±0.01 




-4.2 


-3.8 


0.4720 


23.6 


22.5 


-17.9 


0.4716 


0.81 ±0.01 




-5.5 


7.5 


0.6350 


44.7 


24.1 


-17.1 




<0.02 




-8.9 


3.2 


0.7980 


49.9 


22.4 


-19.5 


0.7969 


0.45 ±0.01 




-3.3 


9.1 


0.5650 


44.1 


23.8 


-17.1 




<0.01 




-6.7 


9.7 


0.3180 


38.3 


20.2 


-18.9 


0.3172 


0.50 ±0.01 




1.6 


14.2 


0.6560 


69.6 


22.7 


-18.6 


0.6556 


1.43 ±0.01 




-12.3 


9.9 


0.7020 


79.1 


21.9 


-19.6 


0.7023 


0.04 ±0.01 




7.9 


-17.0 


0.8280 


99.7 


24.2 


-18.3 




<0.01 




-21.3 


-6.2 


0.3680 


79.2 


19.9 


-19.9 


0.3677 


0.25 ±0.01 




-24.1 


5.3 


0.3680 


81.0 


23.4 


-16.4 








25.5 


17.6 


0.2610 


87.5 


21.6 


-17.3 




<0.03 




12.7 


-31.2 


0.2800 


100.1 


22.9 


-15.9 




< 0.02 



"The R magnitudes presented here were measured in HST/WFPC2 F702W filter for all galaxies but the two in the field around HE0226-41 10. 
No filtered HST images are available for the field of HE0226-41 10. The magnitudes were determined based on a ground-based image obtained 
with the Magellan Baade telescope using IMACS and the Johnson R filter. 



the same dark matter halos of the LRGs. 

To examine the existence of a fiducial relationship between 
the strength of Mg II absorption and galaxy impact parame- 
ter, we present in Figure 3 the distribution of W(2796) with 
respect to p. Circles represent elliptical or SO galaxies, tri- 
angles represent early-type spiral galaxies, and squares repre- 
sent late-type spiral galaxies. Galaxy morphological types are 
adopted from the classifications of Chen et al. (1998, 2001a) 
that were based on a 2D surface profile analysis of the high 
resolution HSTAVFPC2 images. We have also classified the 
galaxy at (-7.9", -7.9") and the galaxy at (11. 6", -12.5") 
from the QSO HE 0226-41 10 as a late-type and an early- 
type spiral galaxy based on their respective spectral features. 
Points with arrows in Figure 3 indicate 3er upper limits. 

Figure 3 shows that there is a moderate anti-correlation be- 
tween W(2796) and p, i.e. weaker Mgll absorption strength 
is observed for galaxies at larger impact parameter. We note 
that the obvious outlier of a Mgll absorber with W(2796) = 
1.43 ±0.01 A at z = 0.656, currently assigned to a galaxy 
at p = 69.6 h~ l kpc, is a known damped Lya absorber of 
log 7V(HI) = 20.36 ±0.08 (Steidel et al. 1997; Rao et al. 
2006). Although the galaxy satisfies our sample selection 
criterion, it is unlikely that the galaxy at p = 69.6 h~ l kpc 
is physically related to the absorber due to the large neutral 
gas content for two reasons. First, all radio detected sources 
in non-cluster environment from deep 21cm surveys have co- 
incident optical counterparts (e.g. Doyle et al. 2005). Sec- 
ond, deep 21cm observations of the M31 halo reveal no HI 
clouds of N(UT) > 10 20 cm" 2 (e.g. Thilkeret al. 2004). This is 
also roughly consistent with the expected contamination rate 
at this projected separation presented in Figure 1. 
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FIG. 3. — Comparison of the corresponding rest-frame absorption equiv- 
alent width W(2796) versus galaxy impact parameter p for a sample of 23 
galaxies. Circles represent elliptical or SO galaxies, triangles represent early- 
type spiral galaxies, and squares represent late-type spiral galaxies. Points 
with arrows indicate 3<r upper limits. 
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To assess the significance of this anti-correlation, we per- 
form a generalized Kendall test that accounts for the pres- 
ence of non-detections and find that the null hypothesis, in 
which W(2796) is randomly distributed with respect to p, can 
be ruled out at a > 97% confidence level. However, the scatter 
of the W(2796) vs. p anti-correlation is clearly large. Specifi- 
cally, 10 of the 23 galaxies at p < 100 h kpc do not produce 
a corresponding Mg II absorber to a sensitive upper limit, im- 
plying a mean covering factor of the Mg + ions (k) ps 57 % 
at p < 100 h~ l kpc. At p < 45 h~ l kpc, we find (k) ps 50 %, 
similar to the preliminary result reported in Tripp & Bowen 
(2005). 

We note, however, that the rest-frame B-band magnitudes 
of the galaxies range from M B -5 log h = -15.9 to -20.3 with 
a median of Mb -5 log h = -19.0, corresponding to a range 
in intrinsic fi-band luminosity from 0.03 Lb, to 1.6 Lb, with 
a median of 0.5L B „ for M B , -5 log h = -19.8 at z ga i ~ 0.4 
(e.g. Faber et al. 2007). Because more massive galaxies are 
expected to reside in more massive halos and therefore pos- 
sess larger gaseous halos, we examine whether accounting 
for galaxy intrinsic luminosity can improve upon the anti- 
correlation displayed in Figure 3. 

Following TC08, we adopt an isothermal profile of finite 
extent R gds for representing the density of the Mg + ions around 
individual galaxies. The expected Mg II absorption equivalent 
width, W as a function of galaxy impact parameter can be 
calculated following 



W(2796) = 



W 



IR 2 . 



p 2 + a 2 



tan 



-P 1 



p 2 + al 



(1) 



at p < R gds and W(2796) = otherwise. The core radius a/, 
is defined to be a/, = Q.2R gds and does not affect the expected 
W(2796) at large p. The extent of Mg II absorbing gas scales 
with the luminosity of the absorbing galaxy according to 



R 



gas 



R„ 



L b , 



(2) 



Following the definition of Equation (2), /? gas * characterizes 
the gaseous extent of L* galaxies. The gaseous extent of 
fainter/more luminous galaxies differs from R giiw according 
the scaling relation. 

To determine the values of Wo, R gd! ,*, and (3 that best de- 
scribe the data, we perform a maximum-likelihood analysis 
that includes the upper limits in our galaxy and Mg II absorber 
pair sample. The likelihood function of this analysis is defined 



m r-oc 



yi-y(Pi,L Bi ) 



(3) 



y'-y(Pi,LB t ) 



where y-, = log W,-, y = log W, and er, is the measurement un- 
certainty of yi. The first product of Equation (3) extends over 
the n measurements and the second product extends over the 
m upper limits. (This definition of the likelihood function is 
appropriate if the residuals about the mean relationship are 
normally distributed.) Because er,- may include significant in- 
trinsic scatter (which presumably arises due to intrinsic vari- 
ations between individual galaxies) as well as measurement 
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FIG. 4. — Comparison of the rest-frame absorption equivalent width 
W(2796) versus galaxy impact parameter p scaled by galaxy B-band lumi- 
nosity. Symbols are the same as those in Figure 3. The solid curve represents 
the best-fit model based on an isothermal density profile of the Mg + ions, 
while the dotted curve shows the expectation from an NFW density profile. 
The errorbar in the upper-right corner represents the intrinsic scatter of the 
best-fit model, in comparison to the observations, as defined in Equation (5). 



error, cr,- is taken to be the quadratic sum of the cosmic scatter 
o> and the measurement error a m 



a. -a c +a mj , 



(4) 



where the intrinsic scatter is defined by 
/ 



= med 



I yi-y(Pi,L B ,)-- [yj-y(pj,L B] )] 
{ 



(5) 

Because _a c depends on the maximum-likelihood solution 
y = logW(pi,LB : ), the maximum-likelihood solution is ob- 
tained iteratively with respect to equations (1) and (2). The 
maximum-likelihood analysis yields 



(3 = 0.35 ±0.05, 



and 



^ g as*=91^/!- 1 kpc 



(6) 



(7) 



and yo = log Wo = -0.25 ± 0.05 for the best-fit 7? gas *. The er- 
rors represent the 1-cr uncertainties of the best-fit parameters. 
In addition, we estimate a cosmic scatter of a c = 0.25 dex. 
The best-fit isothermal model is presented in Figure 4 (solid 
curve). 

The results of our maximum likelihood analysis indicate 
that the extent of Mg II gaseous halo scales with galaxy B- 
band luminosity according to (Lb /Lb ,r an d P = 0.35 ±0.05, 
and that typical L* galaxies possess an extended Mg II halo of 
radius R gds w 91 h~ l kpc. It is worth noting that the best-fit 
9 h kpc characterizes the gaseous extent of galax- 
ies, beyond which no Mg + ions are found to the limit of 
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WX2796) = 0.01 A. Adopting the best-fit scaling relation, the 
gaseous extent of 0.1L* galaxies is expected to be w 41 h~ l 
kpc. Figure 4 shows the relation of W r (2796) vs. luminos- 
ity scaled p, in which seven galaxies at > 100 h~ l kpc do 
not have a corresponding Mg II absorption feature to the limit 
of W(2796) = 0.01 Ang. These non-detections provide the 
strongest constraints on R gd „. 

To examine possible bias in R gd „ as a result of the ini- 
tial sample definition that consider only galaxies at p < 
ro = 100 h~ l kpc, we repeat the likelihood analysis adopting 
ro = 80 h~ l kpc and find that the best-fit results remain the 
same. Choosing a smaller ro reduces the sample substantially 
and the parameters become poorly constrained. Furthermore, 
we note that the results are also insensitive to the choice of 
model gaseous profile. We obtain similar results, adopting a 
Navarro-Frenk- White (NFW) profile commonly used for de- 
scribing the profiles of dark matter halos (Navarro et al. 1997). 
The best-fit NFW profile is also presented in Figure 4 (dotted 
curve). Comparisons of observations and different model pro- 
files show that an isothermal model profile of the underlying 
gaseous halo provides a good fit to the observed W(2796) vs. 
p relation. 

In summary, Figure 4 demonstrates that accounting for the 
scaling relation of gaseous extent with galaxy luminosity, the 
scatter in the W(2796) versus p relation is substantially re- 
duced. Furthermore, the W(2796) versus p relation scaled by 
galaxy intrinsic luminosity allows us to estimate a mean cov- 
ering factor of Mg + ions average across galaxies of a wide 
range of luminosity. We find (k) = 80-86 % within R gds , in- 
cluding the uncertainty in /? gas * from Equation (7). 

5. DISCUSSION 

Using the sample of 13 galaxy-Mg II absorber pairs and 
10 galaxies at p < 100 ft -1 kpc that do not give rise to Mgll 
absorption to a sensitive upper limit, we find that galaxies at 
larger impact parameters produce on average weaker Mgll 
absorbers. This anti-correlation is substantially improved, 
when accounting for the intrinsic luminosity of individual 
galaxies. In addition, there exists a distinct boundary at p = 
R gas , beyond which no Mgll absorbers with W(2796) > 0.01 
A are found. A maximum likelihood analysis shows that 
the observations are best described by an isothermal den- 
sity profile of the gas and a scaling relation R gds = 91^| x 
(Lb /Lb, )(°- 35±005 ) h~ x kpc with a mean covering factor of 
(k) = 80-86 % at p < R gds . The best-fit profile applies to 
galaxies of 0.03 L B , to 1.6L Bt and therefore predominantly 
cold-mode halos. For higher mass halos, we expect from 
TC08 that the corresponding absorption strength for a given 
p is reduced, namely a smaller Wo in Equation (1). In this 
section, we discuss the implications drawn from the results of 
our analysis. 

5.1. Comparisons with Previous Studies 

Adopting the best-fit scaling relation and the gaseous pro- 
file described in Equation (1), we derive a corresponding halo 
size of p w 69 h~ l physical kpc for galaxies of M B -5logh = 
-19.8 at W(2796) = 0.3 A. This halo size is consistent with 
the finding of Lanzetta & Bowen (1990) and the expectation 
inferred from known Mg II frequency function (Nestor et al. 
2005; see the descriptions in § 5.2), but is larger than what was 
reported in Steidel (1995). Steidel (1995) reported a charac- 
teristic size of Mgll halo of R* w 35 h~ l kpc for L» (corre- 
sponding to M B , -5 log h = -19.7 for the cosmology adopted 



in this paper) galaxies, but it is not clear at what W(2796) 
limit the halo size was estimated. At p = 35 h~ l kpc, we ex- 
pect W(2794) = 0.9 A from the best-fit model. 

The high gas covering factor derived from our analysis is 
consistent with what is inferred from the sample of Kacprzak 
et al. (2008). Applying the best-fit luminosity scaling relation 
of gaseous extent, we find that four of the 28 galaxies within 
the expected gaseous extent for W r (2796) > 0.3 A absorbers in 
the Kacprzak et al. sample have a corresponding Mg II weaker 
than the 0.3 A threshold. The inferred covering fraction is 
« 86%. At the same time, this high covering fraction dis- 
agrees with the value reported in Tripp & Bowen (2005), the 
only other on-going survey program designed to constrain gas 
covering fraction by searching for corresponding Mg II ab- 
sorbers at the locations of known field galaxies 6 . The authors 
obtained follow-up spectra of background QSOs that are lo- 
cated within p= 10-39 /T 1 physical kpc of 20 known galaxies 
at Zgai = 0.31-0.55. They found that 10 of the 20 galaxies do 
not produce Mg II absorbers to within a 2-er upper limit of 
W(2796) = 0.1 A, and concluded that the incidence of Mg + 
ions is w 50 % at p < 40 h~ l kpc in individual galactic halos. 

Because more luminous galaxies have larger extended 
Mgll halos according to the scaling relation presented in 
Equations (2), (6), and (7), this discrepancy can be under- 
stood if many of the galaxies in Tripp & Bowen are located 
at impact parameters that are larger than their expected R gds . 
Recall that we obtained a similar estimate of (k) « 57 % at 
p < 100 IT 1 kpc (or (k) « 50 % at p < 45 h~ l kpc), prior 
to accounting for the scaling relation between R gds and L B in 
individual galaxies (§ 4). 

In addition, the sample in Tripp & Bowen includes galaxies 
of luminosity ranging from 0.3 L* to 5L*. Those galaxies that 
are more luminous than L» are also expected to reside in more 
massive halos. At z = 0, 5 L* galaxies are expected to reside in 
> 10 13 h~ x M Q halos (Tinker et al. 2007). If the host dark mat- 
ter halos are substantially more massive than 1O 12 O /i _1 M0, 
then we expect from the halo occupation analysis in TC08 
that they do not contribute significantly to strong Mg II ab- 
sorbers. In summary, further inspections of the properties of 
these galaxies and estimates of the gas covering fraction k as 
a function of radius and W(2796) are crucial for resolving this 
discrepancy. 

5.2. The Origin of Mg II Absorbers 

The scaling relation and the constraint of (k) together also 
allow us to test whether the Mg II absorbers trace a significant 
and representative portion of the galaxy population or merely 
a special class such as those galaxies undergoing a starburst 
episode. If extended Mg II halos are a generic feature of field 7 
galaxies at all epochs, then we can derive the expected number 
density of Mg II absorbers by combining known space density 
of galaxies (versus luminosity) and the product of the cross 
section and covering factor of the Mg + ions. An agreement 
with observational results from surveys of Mgll absorbers 

6 We note two earlier surveys by Bechtold & Ellingson (1992) and Bowen 
et al. (1995) that include a dominant fraction of galaxies in dense cluster en- 
vironment. Both studies reported a substantially lower covering fraction, 
(k) < 20 % for W(2796) > 0.1 absorbers. The low covering fraction of 
strong absorbers is understood, because galaxies clusters reside in massive 
dark matter halos of Mi, > 10 14 /r'MQ and the fraction of "cold" gas in 
these massive halos is expected to be reduced (see TC08). 

7 We define field galaxies as those that are not associated with a known 
cluster. 
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would challenge the scenario that attributes a large fraction 
of these absorbers to starburst systems 8 . 

The predicted number density of Mg II absorption systems 
arising in the extended gaseous halos of galaxies may be eval- 
uated according to 



dN(W > W Hm ) c 



(1+z) 2 



dz 



Ho ^/n M {i+z) 3 +n A 

Lb_ 



(8) 



$(L b ,z)ct(L b ,W;,„,)k, 



where c is the speed of light, <p(Lb,z) is the galaxy luminos- 
ity function, a is the gas cross section for producing Mg II 
absorbers of W(2796) > W;, m that scales with galaxy B-band 
luminosity, and k is the halo covering factor. Substituting the 
scaling relationship according to Equations (2), (6), and (7), 
we find 



dN(W > W Um ) c it (k) [R'{W lim )} 2 (\ +zf 



dz 



Ho ^/Q M (l+z) 3 + n A 



(9) 



Lb, 



Lb_ 
Lb, 



2 6 



®(Lb,z), 



where R' gas *(Wiim) is the extent of Mg II halos that corresponds 
toW(2796) = WW 

To evaluate Equation (9), we adopt the luminosity function 
of the blue galaxy population at z ga ] = 0.3 - 0.5 from Faber 
et al. (2007), which is characterized by Mb, — 5 log = -19.8, 
4>* = 9.6 x 10~ 3 /z 3 Mpc" 3 , and a faint-end slope a = -1.3. 
Including galaxies of Lb > 0.01 Lb, and our best-fit /? gas * and 
(k), Equation (9) yields dN/dz = 0.028-0.03 for W(2796) > 
2 A absorbers, dN/dz = 0.21-0.23 for W(2796) > 1 A ab- 
sorbers, and dN/dz = 0.88-1.07 for W(2796) > 0.3 A ab- 
sorbers at z = 0.5. The predictions based on galaxies more 
luminous than 0.01 L* agree well with the observed dN/dz 
from Nestor et al. (2005), lending strong support for the hy- 
pothesis that the Mg II absorbers trace a significant and repre- 
sentative portion, rather than a sub-class, of the galaxy popu- 
lation. 

The model profile of Equation (1) has been proved to be 
a good representation of the galaxy and absorber pairs at p = 
16.3- 100 h kpc, when the scaling relation of Equations (2), 
(6), and (7) is accounted for. An important test is to examine 
whether the model remains a good fit at smaller p and strong 
absorbers. Identifying galaxy-absorber pairs at small impact 
parameters is difficult, because the glare of the background 
QSO often prohibits identification of faint galaxies at angular 
distances < 2" that corresponds to ~ 10 h~ l physical kpc at 
z<0.5. 

Recent studies of intervening absorption-line systems along 
the lines of sight toward the optical afterglows of distant 7-ray 
bursts (GRB) have uncovered a number of strong Mgll ab- 
sorbers (e.g. Prochter et al. 2006b). Unlike QSOs, GRB after- 
glows disappear after a while, permitting exhaustive searches 
of faint galaxies at 9 < 2" that give rise to strong Mg II ab- 
sorbers along the lines of sight. A redshift survey of galax- 
ies along the sightline toward GRB 060418 (z G rb = 1.491) 

8 However, we will not to able to rule out a scenario in which the blow-out 
winds are common in most galaxy histories. Given enough time, long-lived 
outflowing gas gets mixed in with existing halo gas and newly accreted ma- 
terial from the IGM. There is little distinction at this point between outflows 
and accreted clouds. It is also not clear whether the outflow material would 
maintain the same degree of complex gas kinematics as those that are closer 
to the parent starburst regions. 
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FIG. 5. — The W(2796) versus p correlation scaled by galaxy B-band lumi- 
nosity presented in Figure 4, but including pairs from Pollack et al. (2008; star 
points) and Bouche et al. (2007; points with horizontal errorbars) for compar- 
isons. The solid curves represent the best-fit model based on an isothermal 
density profile, with a/, =0.1 i? gas (top curve) and a/, = 0.2R ga s (bottom curve). 
The pairs of Pollack et al. are for galaxies and Mg II absorbers identified at 
z = 0.6— 1 along the sightline toward GRB 060418. The errorbars in these data 
points are smaller than the size of the symbols. The pairs of Bouche et al. 
are galaxies found at the positions of known Mg II absorbers based on the 
presence of Ha emission lines. Errorbars represent the uncertainties of the 
inferred Mb from the Mb versus L(Ha) correlation published in Tresse et al. 
(2002). Crosses represent the absorbers in the Bouche sample that are likely 
(> 40 % probability) to be a damped Lya absorber of N(H1) > 2 X 10 20 cirT 2 
based on the presence of strong Fe II and Mg I transitions (Rao et al. 2006). 



has uncovered the absorbing galaxies for three strong Mg II 
absorbers at z = 0.603-1.107 (Pollack et al. 2008, in prepa- 
ration), similar to the redshift range covered by our sample. 
The rest-frame Z?-band magnitudes of the galaxies range from 
M B -5 log h = -16.8 to M B -5 log h = -18.3; Impact parame- 
ters of the galaxies range from p = 7.5 h~ l kpc to p = 16.5 h~ l 
kpc. Adopting the best-fit scaling relation, we include the 
three galaxy-Mgll absorber pairs in the W(2796) versus p 
plot for comparison (star points in Figure 5). The good 
agreement between the best-fit model and the absorbers found 
along GRB sightlines provides further support for our model 
at small impact parameters. 

5.3. Are Strong /W(2796) > I A] Absorbers at z ~ 1 
Produced Primarily in Starburst Outflows? 

As described in § 1, the large line width of W(2796) > 1 A 
Mg II absorbers has recently been interpreted as due to out- 
flow motions in starburst driven winds. This scenario was ini- 
tially motivated by the observed inverse correlation between 
W(2796) and the clustering amplitude of Mgll absorbers at 
(z) = 0.6 (Bouche et al. 2006). In TC08, we have shown that 
this inverse correlation can be understood as due to an ele- 
vated clustering amplitude of W(2796) < 1 A absorbers from 
contributions of survived cold clouds in massive hot halos, 
rather than a suppressed clustering strength of W(2796) > 2 
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A absorbers. In addition, we have shown in §§ 4, 5.1, and 
5.2 that extended Mgll halos are a generic feature of field 
galaxies over a wide luminosity range. The predicted number 
density of the absorbers from adopting the scaling relation, 
Equations (2), (6), and (7), and the galaxy luminosity func- 
tion agrees well with observations. Specifically, we expect 
dAf/dz ~ 0.03 for W(2796) > 2 A absorbers, in comparison 
to the observed dAf/dz = 0.03-0.04 for W(2796) > 2 A ab- 
sorbers at z ~ 0.5 from Nestor et al. (2005). While there is no 
clear indicator one can apply to unambiguously distinguish 
between gas inflows and outflows in the distant universe, we 
discuss several caveats related to the starburst outflow sce- 
nario. 

If the absorbers are produced in starburst driven outflows, 
then a large fraction of the absorbing galaxy population might 
be expected to exhibit disturbed morphology (Mobasher et al. 
2004). The Mgll absorbing galaxies in our sample exhibit 
regular disk morphologies in the HST images shown in Fig- 
ure 2. In particular, the strong absorber found with W(2796) = 
1 .55 A at z = 0.892 toward 3C336 is associated with an edge- 
on disk galaxy at p = 16.3 h~ l kpc. The three Mg II absorbers 
of W(2796) > 1 A found toward GRB 0604 1 8 also exhibit reg- 
ular disk structures in high-resolution images obtained both in 
space using HST and on the ground using LGAO on the Keck 
telescopes (Pollack et al. 2008, in preparation). Kacprzak 
et al. (2007) also showed that there is little correlation be- 
tween W r (2796) and galaxy asymmetry for Mg II absorbers of 
W r (2796) > 1.4 A. While some strong Mg II absorbers have 
been presented to show evidence of arising in expanding su- 
perbubble shells (Bond et al. 2001), available high-resolution 
images of other strong Mg II absorbers do not support out- 
flow being principally responsible for the observed strong ab- 
sorbers. This conclusion is also consistent with the spectral 
properties of Mg II-selected galaxies published by Guillemin 
& Bergeron (1997) 

On the other hand, Bouche et al. (2007) presented new ob- 
servations from a targeted search for corresponding Ha emis- 
sion in the vicinity of known Mg II absorbers. These authors 
successfully identified Ha emission near 14 of 21 Mgll ab- 
sorbers with W(2796) > 1.9 A at z ~ 1. The integrated raw 9 
Ha luminosities range from L^[ a = 6.8 x 10 40 h erg s" 1 to 
Ljjq, = 8.4 x 10 41 h~ 2 erg s" 1 . The impact parameters range 
from p = 1.4 h~ x kpc to p = 35 h~ l kpc. Adopting a mean 
extinction correction A v = 0.8 mag, the authors inferred a 
mean star formation rate of (SFR) = 5 .9 h~ 2 M Q yr" 1 for these 
galaxies. 

The Ha luminosity function at z ~ 1 has been studied 
by a number of groups (see Doherty et al. 2006 for a list 
of references). It is characterized by a Schechter function 
of L Ha = 5 x 10 41 h 2 erg s" 1 , = 0.013 h 3 Mpc~ 3 , and 
a = —1.3. The galaxies detected in Ha by Bouche et al. cor- 
respond to 0. 14 - 1 .68 Lj| a at z ~ 1 , and therefore have only 
modest L(Ha). Integrating the Ha luminosity function, we 

9 No extinction correction is applied to L(Ha) in our analysis, because the 
dust content is expected to vary substantially in z ~ 1 galaxies (e.g. Tresse 
et al. 2002) and is unknown for galaxies in the Bouche et al. sample. Bouche 
et al. (2007) applied a constant correction to all observed Ha flux in their 
calculation, resulting in a simple scaling of the observed values. The scatter 
in their measurements induced by dust remains the same. For consistency, 
we consider in the following discussion only observed values without extinc- 
tion correction both for the absorbing galaxies and for comparison results, 
such as the Ha luminosity function and the Mg-L(Ha) correlation, from the 
literature. 



obtain a mean space density of 0.034 h 3 Mpc 3 for galaxies 
of brighter than 0.1 Ljj a . This is already 50% of all galaxies 
more luminous than 0.01 L B , at z ~ 1, according to the best- 
fit luminosity function of the blue galaxy population from 
Faber et al. (2007). This exercise shows that Ha emission is 
commonly seen in z ~ 1 galaxies and the presence of modest 
L(Ha) does not argue for the presence of starburst outflows, 
unless the outflow material is long-lived. 

For comparison, we include the 14 galaxy and absorber 
pairs from Bouche et al. (2007) in the W(2796) versus p plot 
in Figure 5, applying the same scaling relation described in 
Equations (2), (6), and (7). We estimate Mb of these galaxies 
based on the published M B -L(Ha) correlation in Tresse et al. 
(2002). We adopt the observed L(Ha) and the scaling relation 
Mb = 73.5-2.27 x logL(Ha) that best describes the data in 
Tresse et al. (2002). At small impact parameters, the model 
predictions depend sensitively on the choice of a/, in Equation 
(1). In addition to the fiducial model, we also include in the 
figure a model with a h = 0.1 /? gas . 

Nine of the 14 Mg II absorbers are considered likely (> 40 
%) damped Lya absorbers (DLAs) of JV(HI) > 2 x 10 20 cm -2 
based on the presence of strong Fell and Mgl transitions 
(Rao et al. 2006). DLAs are the high-redshift analog of neu- 
tral gas regions that resembles the disks of nearby luminous 
galaxies (e.g. Wolfe et al. 2005). We therefore expect that 
the impact parameter distribution of galaxy-DLA pairs rep- 
resents the typical extent of H I disks at high redshifts. We 
highlight the potential DLAs in the Bouche et al. sample as 
crosses in Figure 5. After applying the best-fit scaling rela- 
tion for Mg II absorbers, we find that the impact parameters 
of these galaxy and potential DLA pairs range from 8 h~ x kpc 
to as high as 38 h kpc. The three largest separation pairs 
have p' > 28 h~ x kpc, where p' is the luminosity scaled pro- 
jected distance. The scaling of impact parameter is justified, 
because smaller/fainter galaxies have on average smaller HI 
disks (e.g. Cayatte et al. 1994). The large impact parameters, 
together with the likely high A^(HI), suggest that the galaxies 
responsible for these absorbers (presumably at smaller angu- 
lar distances from the QSOs) may still be missing. For the re- 
maining sample, our model is in reasonably good agreement 
with their data. 

In summary, we find a lack of empirical measurements that 
argue for starburst outflows to be a dominant mechanism for 
producing the observed strong Mg II absorbers at z ~ 1 . On 
the other hand, both the observed number density of Mg II ab- 
sorbers and the W(2796) versus p correlation can be explained 
by extended Mg II halos around typical field galaxies. The 
fraction of Mg II absorbers originating in starburst outflows is 
therefore expected to be small at z ~ 1, unless starburst out- 
flows are a common feature in field galaxies. We caution, 
however, that the same conclusion cannot be automatically 
applied to absorbers at z ~ 2 before a similar analysis is per- 
formed. 

5.4. Constraints on the Extended Gaseous Halos Around 
Galaxies 

The sharp decline of W(2976) found at p w 90 h~ x kpc in 
Figure 4 clearly indicates that Mg + ions have a finite extent in 
the extended gaseous halos around galaxies. This is reminis- 
cent of the sharp boundary between the C IV absorbing and 
non-absorbing regions around (z) = 0.4 galaxies reported in 
Chen et al. (2001b). 

Given the nature of a photo-ionized gas (Bergeron & Stasfn- 
ska 1986; Hamann 1997), we examine whether the finite ex- 
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FIG. 6. — Ionization fractions of Mg + ions and H° atoms versus total gas 
density «h for photo-ionized gas. The values are calculated using the Cloudy 
software (Ferland et al. 1998; version 06.02) for gaseous clouds of plane par- 
allel geometry and 100 pc thickness. The input ionizing radiation intensity is 
set to be /„(912) = 7 X 10~ 23 erg s" 1 Hz" 1 cm" 2 Sr 1 suitable for the z < 1 
universe (Scott et al. 2002). We perform the calculations, assuming solar 
(dotted curve for Mg and dot-dashed curve for H) and 1 / 10 solar metallicity 
(solid curve for Mg and dashed curve for H). 



tent of Mg + ions can be understood as due to photo-ionization 
of the absorbing clouds at large radii. In a two-phase halo 
model (e.g. Mo & Miralda-Escude 1996), cold clouds that are 
responsible for producing the absorption features are pressure 
confined in a hot medium. As the total gas density declines to- 
ward larger radii, the density of the cold clouds also decreases 
and the clouds become optically thin to the ultraviolet photons 
that can quickly ionize the Mg + ions (LP. = 15.035 eV). With 
a known background radiation field, the location where the 
ionization transition occurs may be adopted to constrain the 
density of cold gas at large radii. 

We carry out the photo-ionization analysis, using the 
Cloudy software (Ferland et al. 1998; version 06.02), and cal- 
culate the ionization fractions of Mg + ions for a grid of models 
with different total gas density and metallicity. We adopt the 
spectral shape of the UV background radiation field of Haardt 
& Madau in Cloudy for the incident ionizing flux, and scale 
the radiation intensity to the background UV radiation field 
inferred from QSO proximity effect (Scott et al. 2002). We 
adopt a plane parallel geometry for the clouds with a thick- 
ness of 100 pc, which is motivated by the observed coherent 
length of Mg II absorbers toward lensed QSOs (Rauch et al. 
2002). The slab of gas is illuminated by the ionization radia- 
tion field on both sides. 

Figure 6 shows the expected fractions of Mg + ions (solid 
and dotted curves) as a function of the total gas density «h- 
We include the calculated fraction of H° atoms (dashed and 
dash-dotted curves) for comparison. The solid and dashed 
curves are for gas of 1/10 solar metallicity, while the dotted 
and dash-dotted curves are for gas of solar metallicity. Our 
photo-ionization calculations indicate that indeed the fraction 
of Mg + declines linearly at n H ^ 0.025 cm" 3 for a mean radi- 



ation intensity of J v {9\2) = 7 x 10" 23 erg s" 1 Hz" 1 cm" 2 Sr" 1 
(which is appropriate for the observed background radiation 
field at z < 1 Scott et al. 2002). However, the transition is 
not as sharp as what is seen for hydrogen atoms. Including 
the estimated ionization fraction from the Cloudy models and 
adopting the density profile parameterized in Equation (1), 
we derive the expected correlation between W r (2796) and p 
and find that photo-ionization is insufficient for explaining the 
presence of the sharp decline in the observed Mg II absorption 
strength. 

This extent of Mg II gaseous halos (« 90 h~ x kpc around 
typical L* galaxies) is remarkably similar to the extent of 
CIV gaseous halo (« 110 h~ l kpc around typical L* galax- 
ies, after correcting for the A-cosmology) reported by Chen 
et al. (2001b). On the basis of 14 galaxy and CIV absorber 
pairs and 36 galaxies that do not produce corresponding C IV 
absorption lines to within sensitive upper limits, these au- 
thors find that the extent of C IV-absorbing gas around galax- 
ies scales with galaxy B-band luminosity as R gas (CTV) oc 
110 x L^' 5±01 /T 1 kpc (corrected for the A cosmology). In 
addition, there exists a sharp boundary between C IV absorb- 
ing and non-absorbing regions at /? gas (CIV). But the scatter in 
the observed CIV absorption strengths at p < /? gas (CIV) be- 
tween different galaxies appears to be significantly larger than 
what is seen in Figure 4, and there is a lack of anti-correlation 
between the observed C VI absorbing strength and galaxy im- 
pact parameter. 

Given the distinct ionization potentials of the C 3+ and Mg + 
ions, the location where the abrupt transition between absorb- 
ing and non-absorbing regions occurs can be understood as 
a critical radius below which cool clouds can form and sta- 
blize in an otherwise hot halo. This two-phase model to in- 
terpret QSO absorption line systems was formulated in Mo & 
Miralda-Escude (1996) and later re-visited by Mailer & Bul- 
lock (2004). 

Taking into account the appropriate photo-ionization condi- 
tion of the cold clouds, Mo & Miralda-Escude presented the 
expected Mg II and C IV absorption strength versus radius for 
halos of different mass. Their models showed a lack of radial 
dependence on the C IV absorption strength and a tight corre- 
lation between Mg II W r (2796) and p for halos of a wide range 
of mass. The differences in the radial dependence of different 
ions is understood as being due to the photo-ionization con- 
dition of cold clouds in the UV background radiation field. 
These model expectations agree with observational findings 
very well. 

The agreement between observations and a simple two- 
phase model is encouraging, although it is not clear how the 
two-phase model applies to cold-mode halos (e.g. Dekel & 
Birnboim 2006). Measurements of the relative abundances 
between C 3+ and Mg + as a function of radius will provide ad- 
ditional support for the origin of these metal-line absorbers, 
further constraining the metallicity and density of halo gas 
around galaxies. 

5.5. Implications for the Halo Occupation Analysis and 
Future Work 

Our initial halo occupation analysis (TC08) included two 
principal assumptions for the extent of Mg + . First, we as- 
sumed that the extent of the cold gaseous halo is /? gas ,i2 = 
50/T 1 physical kpc for M h = l0 n h~ 1 M Q halos, which is 
roughly 1 /3 of the halo radius /?200 at z = 0.6. Recall that 
the standard definition of a dark matter halo is an object with 
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a mean interior density of 200 times the background density. 
Second, the gaseous extent scaled with the halo mass as M' h 
and t = 1/3, following the expectation of the scaling relation 
between halo mass and the corresponding virial radius. The 
results presented in § 4 offer a direct test of these assumptions. 

To compare the empirical constraints of Equations (6) & 
(7) with our initial model assumptions, we first derive the cor- 
responding dark matter halo mass M/, and their halo radius 
/?200 for galaxies of known absolute B-band magnitudes. The 
adopted M B « for the analysis in § 4 corresponds to a mag- 
nitude of Mb, -5 log h = -19.7 in the bandpass of the Two- 
Degree Field Galaxy Redshift Survey (2dFGRS, Colless et al. 
2001). The 2dFGRS probes galaxies at z ~ 0.1. Excluding 
galaxies that exist as satellites in a cluster environment 10 , the 
mean halo mass forM B „ galaxies is lO 125 /r 1 /r 1 M at z ~ 
0.1 (Tinker et al. 2007; van den Bosch et al. 2007). At z ~ 1, 
galaxies of this magnitude on average reside in 10 n - 9 /i _1 M Q 
halos (ZCZ07). Interpolating in log (1 +z), we find that galax- 
ies of M B , -5 log h = -19.8 are expected to reside in halos of 
M h = \0 l23 h~ l M Q at z = 0.4 and that the halo radius forM B „ 
galaxies is R 200 = 212 h~ l physical kpc. The best-fit /? gas * from 
Equation (7) therefore implies /? gas w 0.4/? 200- 

Next, we derive the expected scaling relation between /? gas 
and halo mass M h . The halo occupation analysis of 2dFGRS 
galaxies from Tinker et al. (2007) showed that the relation- 
ship between Mh and luminosity for galaxies of Lh, < L* at 
Z ~ 0.1 is M h = 10 12 ' 5 (L B /L B J 1 3 h~ l M . this monotonic re- 
lationship is appropriate for galaxies that reside at the centers 
of their dark matter halos and are the brightest galaxy in the 
halo. A similar scaling relation was obtained from DEEP2 
data by ZCZ07 for galaxies at z ~ 1 . Because there is a mono- 
tonic relationship between halo mass and galaxy luminosity, 
the scaling relation of R g . ds with L B in Equation (2) is equiva- 
lent to a scaling relation between Mh and B200 that are related 
by Mh oc ^200- Adopting (3 = 0.35 from Equation (6), we de- 
rive R gd>i = R gd „ x [Mh/(10 l23 h- l M G )] 03 h~ l kpc at 0.4. 

The calculations above show that observations support the 
initial assumptions in TC08 that the gaseous radius is a con- 
stant fraction of the halo radius for all halos and that the 
gaseous radius scales with halo mass according to M' h with 
t 1/3. We note that our galaxy sample therefore probes halo 
mass range from 10 10 ' 5 /z _1 M Q to 10 12 ' 8 /;" 1 M Q , and that the 
best-fit scaling relation is driven by galaxies in the cold-hot 
transition regime. Next, we compare the observed covering 
fraction of Mg II absorbers with the predicted mean of TC08. 

As summarized in § 1, we find in the halo occupation anal- 
ysis of TC08 that dark matter halos of M h = 10 11 5-12 5 h~ l M 
on average have a covering fraction of unity (k(M/,)) = 1 
at R gas < /?20o/3 for Mg + and (n{M h )) = at larger radii. 
This halo mass range corresponds to a luminosity range of 
L B = 0.2 - 1 ALb, at z ~ 0.4, following the mass-to-light scal- 
ing relation discussed above. Sixteen galaxies in our sample 
have L B > 0.2L B ,, eleven of which are at log p+0.14 x (M B - 

10 We note that the fraction of satellite galaxies is < 15 % for L» galaxies 
(Tinker et al. 2007; ZCZ07). In addition, none of the galaxies in our sample 
are observed to be part of a larger group or cluster. Here we carry out our 



Mb ,) < log /? gas *. We identify a corresponding Mg II absorber 
for every one of these galaxies, indicating a mean covering 
fraction of (k) = 100 % at R < R g . d>i . This empirical result 
agrees with our model expectation. To assess the significance 
of detecting a Mg II absorber in every galaxy in the sample of 
eleven, we adopt the binomial distribution function and find 
that we can constrain the underlying mean covering factor at 
(k) > 0.76 with greater than 95% confidence level. 

Three of the seven remaining L B < 0.2L B , galaxies in our 
sample have log p + 0.14 x (M B -M B J < log/? gas ». Two of 
them do not have corresponding Mg II absorbers identified to 
a sensitive upper limit, implying (k) = 33 % at R < /? gas for 
fainter galaxies (presumably lower-mass halos). This is also 
consistent with the expectation of TC08, where we found a 
quickly declining (k) for halos of M h < 10 115 h~ l M Q , al- 
though the observations have a much lower statistical signifi- 
cance. In addition, we note that there exists an intrinsic scatter 
in the fiducial scaling relation between Mh and Lb- This intrin- 
sic scatter is expected to smear out the exact edge of gaseous 
halos, particularly when considering an ensemble of objects 
with a wide range of luminosity. 

In summary, the initial results from the on-going survey 
of Mg II absorbers in the vicinity of known galaxies close 
to QSO lines of sight have allowed us to study in detail the 
properties of extended gas around galaxies. The best-fit scal- 
ing relation between the extent of Mg + and galaxy luminos- 
ity and the observed mean gas covering fraction together pro- 
vide strong empirical evidence to support the halo occupa- 
tion model of TC08. To improve the statistical significances 
of various empirical parameters and to obtain a better under- 
standing of the incidence of cold gas in dark matter halos, 
such as k vs. Mh and k vs. p, a larger sample of galaxy and 
Mg II absorber pairs established for a wide range of galaxy 
luminosity and morphology is necessary. In particular, only 
three galaxies in our current sample would be considered in a 
"cold-mode" halo based on their luminosity and none shows 
a corresponding Mgll absorber. Deep surveys targeting at 
fainter dwarfs are necessary to probe these "cold-mode" ha- 
los. Finally, we encourage observations of Mg II absorbers 
along the lines of sight toward close QSO pairs for obtaining 
a direct measure of the incidence of Mg + in individual halos. 
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calculations, assuming that most of galaxies in our sample are central galaxies 
in their dark matter halos. 
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